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ABSTRACT 

This  paper  intends  to  provide  a  description  of  the  LAREDO  tool  including  architecture,  design  choices 
and  mathematical  formulations.  The  goal  of  LAREDO,  which  is  a  six  DOF  simulator  developed  under 
Matlab/Simulink-based  environment,  is  to  reproduce  the  Rendezvous  and  Docking  (RvD)  sequences 
applicable  to  different  exploration  mission-scenarios  and  to  provide  detailed  analysis  of  formation  flyings 
(FF)  satellite  constellations  using  highly  coupled  GNC  algorithms,  which  integrate  attitude,  navigation 
and  maneuvers.  Among  those,  it  is  possible  to  select:  Earth,  Moon,  Mars,  Venus,  Deimos,  Triton  and 
Europa  as  final  destination.  Within  the  LAREDO  simulator  a  launching  optimization  module  is 
implemented  to  optimize  the  launching  phase  once  it  has  lifted  off  from  the  Earth,  Moon  or  Mars  surfaces. 
This  optimization  is  performed  such  that  maximum  payload  is  delivered  into  orbit.  The  tool  simulates  the 
rocket  stages  and  also  evaluates  the  launch  opportunities  for  optimum  launch.  The  orbit  control  can  be 
performed  by  impulsive  or  continuous  maneuvers  both  on  circular  and  elliptical  orbits.  The  type  of 
maneuvers,  their  order  and  schedule  is  completely  configurable  by  the  users.  The  maneuvers  are  very 
accurate  even  under  perturbations.  In  case  the  identified  safety  limits  are  trespassed,  a  collision 
avoidance  maneuver  (CAM)  will  be  automatically  activated,  moving  the  spacecraft  to  a  safety  region.  The 
simulator  propagates  considering  various  orbit  perturbations  (aerodynamic,  non-sphericity,  solar 
radiation  and  third  body  perturbations),  communications  (both  with  ground  stations  and  between 
satellites),  and  provides  full  analysis  of  the  trajectory.  The  GNC  modules  incorporate  sensor  models  for 
LIDAR,  camera,  IMU  and  RF  devices,  and  both  orbit  and  attitude  guidance  and  control.  Furthermore, 
LAREDO  offers  powerful  3-D  visualization  of  trajectories  and  attitude,  as  well  as  other  important  features 
such  as  consumed  A  V,  ground  station  contacts,  eclipse  times,  and  others. 


1.0  NOMENCLATURE 


The  following  concepts  and  representations  are  used  across  this  paper: 


PAV 

CRF 


Chaser 

Target 


=  Planetary  Ascent  Vehicle.  In  the  launch  ascent  context,  this  is  the  spacecraft  that  takes  off 
from  the  planetary  surface  and  is  inserted  into  orbit. 

=  Cylindrical  Reference  Frame.  Y-axis  is  Perpendicular  to  the  orbital  frame,  Z-axis  is 
opposite  to  the  Target  position  vector  in  the  MEE2000  inertial  frame,  and  the  X-axis  is 
orthogonal  to  the  Z  and  Y. 

=  In  the  context  of  satellite  formations,  the  Chaser  is  any  of  the  satellites  in  the  formation 
that  is  not  the  Target  one. 

=  In  the  context  of  satellite  formations,  guidance  and  control  schemes  use  the  Target  satellite 
motion  as  the  reference,  fixing  with  respect  to  it  the  nominal  positions  and  velocities  of 
the  Chaser  satellites. 
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Launch  ascent  optimization  variables 
Payload  mass 

Ignition  time  of  the  jth  rocket  stage 

Actual  simulation  time 

Initial  and  final  position  vectors 

Initial  and  final  reference  position  vectors 

Initial  and  final  velocity  vectors 

Initial  and  final  reference  velocity  vectors 


to*1/ 

Av0 ,  Av1 
xc,xt 

CO 


Initial  and  final  orbital  maneuvers  times 
Initial  and  final  orbital  maneuver  impulses 
Chaser  and  Target  spacecraft  state  vectors 
Orbital  angular  rate. 

Acceleration  in  the  x,  y  and  z-axis. 


2.0  ARCHITECTURE 

The  architecture  of  the  software  is  represented  in  the  Fig.l.  In  this  figure  it  can  be  observed  that  the 
Launcher  Optimization  and  the  RvD  parts  run  separately. 

However,  it  has  to  be  stressed  that  the  RvD  part  can  import  the  last  state  of  the  Launcher  Optimization  part 
to  start  a  RvD  simulation.  The  simulator  contemplates  orbital  and  attitude  Guidance,  Navigation,  and 
Control,  as  well  as  sensors  and  actuators.  The  outputs  from  both  the  Launcher  Optimization  and  RvD  parts 
can  be  visualized  in  the  Figures  Of  Merit  (FOM)  section. 


Figure  1  General  Simulator  Architecture. 


3.0  LAUNCHER  OPTIMIZATION 

LAREDO  tool  is  able  to  simulate  the  ascent  trajectory  of  a  vehicle  from  the  surfaces  of  the  Earth,  the 
Moon  and  Mars.  The  software  is  able  not  only  to  simulate  the  dynamics  of  the  rocket,  but  also  to  optimize 
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the  trajectory  to  reach  the  desired  final  orbit  with  the  required  relative  argument  of  latitude.  The  trajectory 
propagation  considers  the  atmospheric  and  the  non-sphericity  perturbations,  as  well  as  the  rocket  stages 
and  the  fuel  mass  losses. 

The  rocket  has  two  configurable  stages,  where  the  structural  masses,  propellant  masses,  and  other 
parameters  can  be  defined.  A  coast  arc  can  be  inserted  between  the  two  stages,  where  the  duration  is  an 
optimization  variable. 

3.1  Rocket  dynamics  modeling 

One  of  the  key  parameters  when  modeling  a  rocket  ascent  trajectory  is  the  rocket  thrust  vector.  The  thrust 
vector  is  defined  by  a  modulus  and  a  direction,  given  by  Eq.(l)1  and  Eq.(2)  respectively: 

dm^  j  j 

=  '  SO  I sp  ~  (1) 

9  is  the  sea  level  gravity,  /  J  is  the  specific  impulse  in  vacuum 
area  of  the  jth  stage  nozzle  and  Pa  is  the  atmospheric  pressure. 

j  AJ'  +  BJ  [tj  —fj  start  ) 

7  -  7  .  - : - (  (2) 

AJ  +  B  ^  \ti  —t^ start  J| 

Eq.(2)  which  expresses  the  thrust  direction,  is  also  known  as  the  44  bi-linear  tangent  law2  ”.  Bothv42an(j 
BJ  vectors  are  optimized  when  the  second  stage  of  the  rocket  is  burning.  During  the  first  stage,  the  rocket 
direction  is  considered  to  remain  constant  and  thus  BJ  =  0. 

The  mass  of  the  burning  stage  is  continuously  reduced.  The  loss  of  mass  will  be  determined  by  the  rocket 
flow  rate  by  the  following  expression: 


\t\j 


dm 


where  is  the  jth  stage  flow  rate,  g{ 
dt 

of  the  jth  stage,  AtJ  is  the  effective  exit 


f 


J  _ 


m 


ti-l 


dm  J 
dt 


(3) 


dm  ^ 

where  -  is  the  jth  stage  flow  rate,  mti_/  is  the  jth  stage  mass  at  the  previous  simulation  step  time  A/ 

dt 


dm 


and  mti  is  the  j  stage  mass  at  the  time  tt.  Note  that -  will  be  zero  if  the  stage  is  not  burning. 

dt 


If  4V’  is  considered  as  the  total  number  of  stages,  the  mass  of  the  rocket  is  defined  as: 


mi 


-  ZV  +  m  payload 

j=l 


(4) 
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The  heat  transfer  rate  is  one  of  the  main  factors  to  be  considered  when  launching  a  vehicle.  If  the  heat 
transfer  exceeds  certain  maximum  allowable  value,  the  rocket  and/or  its  cargo  and  crew  can  suffer  severe 
damage  or  even  disintegrate.  The  heat  transfer  equation  is  given  by  the  following  equation  (Sutton  & 
Graves  correlation3): 


q  =  C1  x 


(  \C2 

Pa 

\^n  J 


un 


1000 


(5) 


where  p ^  is  the  atmospheric  density,  U ^  is  the  velocity  in  the  corotating  reference  frame  (reference 

frame  fix  to  the  planet),  Rn  is  the  rocket  nose  radius  and  Q,  C2  and  C3  are  coefficients.  For  the  case  of 
Mars,  the  C;,  C2  and  C3  coefficients  are  18.9,  0.5  and  3  respectively. 

3.2  Optimization 

The  goal  of  most  rockets  is  to  achieve  the  required  orbit  with  the  maximum  injected  payload  mass.  The 
software  tries  therefore,  to  minimize  the  following  cost  function: 

J  =  ~m  payload  (6) 

The  trajectory  of  the  rocket  is  propagated  using  a  full  Cowell  propagator  with  atmospheric  and  non¬ 
sphericity  perturbations  acting  upon  the  rocket,  and  also  taking  into  account  the  different  rocket  phases. 

The  cost  function  expressed  in  Eq.(6)  is  minimized  while  satisfying  a  few  equality  and  inequality 
constraints  which  ensure  that  the  trajectory  is  feasible. 

•  Equality  constraints:  The  equality  constraints  are  defined  by  the  desired  final  orbital  parameters. 
This  way,  the  orbital  parameters  propagated  for  the  launch  ascent  with  the  Cowell  full-  propagator 
must  be  equal  to  the  final  desired  orbit  parameters. 

•  Inequality  constraints:  Three  inequality  constraints  are  imposed.  The  first  one  is  to  make  sure  that 
the  height  achieved  when  the  launch  optimization  is  performed  is  greater  than  the  planet  radius. 
The  second  one  is  to  ensure  that  the  heat  rate  will  never  exceed  a  maximum  value.  Finally,  the 
third  one  is  to  ensure  that  the  PAV  spacecraft  is  placed  with  respect  to  the  other  satellite  within  a 
certain  relative  argument  of  latitude  limits. 

4.0  ORBITAL  MANOEUVRES  FOR  RVD 

The  use  of  more  than  one  spacecraft  to  perform  a  mission  makes  the  RvD  maneuvers  more  and  more 
important.  These  maneuvers  are  usually  performed  in  circular  orbits,  however,  there  is  rising  interest  in 
RvD  in  elliptical  orbits.  LAREDO  is  able  to  provide  a  very  accurate  analysis  in  both  circular  and  elliptical 
orbits,  estimating  the  AV,  the  fuel  consumption  and  other  important  features. 

RvD  is  a  very  risky  part  of  the  mission  as  any  failure  can  cause  the  satellites  damage.  The  tool  simulates  a 
Collision  Avoidance  Maneuvers  that  is  activated  if  the  satellites  trespass  certain  safety  limits.  These  limits 
are  configurable,  and  of  course  dependent  on  the  relative  satellite  distances.  When  the  satellites  are  far 
away  from  each  other  the  safety  zone  is  assumed  to  be  a  sphere,  whereas  if  they  are  close  the  safety  zone 
is  modeled  as  a  cone  (approaching  corridor).  If  the  satellite  performs  some  of  the  programmed  maneuvers 
crossing  any  of  these  boundaries,  the  required  orbital  and  attitude  maneuvers  are  automatically  activated. 
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4.1  Corrective  maneuvers  for  circular  orbits 

The  most  important  problem  in  RvD  comes  out  from  the  high  maneuvering  precision  required.  This  severe 
requirement  during  the  RvD  implies  reducing  as  much  as  possible  the  perturbations  effects  acting  upon  the 
spacecraft.  In  order  to  counteract  these  perturbing  effects,  the  spacecraft  will  perform  corrective 
maneuvers  during  the  whole  trajectory.  An  increase  in  the  number  of  corrective  maneuvers  generally 
implies  a  reduction  in  the  final  trajectory  error,  but  with  a  penalty  on  the  consumed  fuel  (see  Fig.  2). 


Vf,ref 

Vf,ref 


v0+Av0 

h 


Figure  2.  Corrective  maneuvers. 


The  method  implemented  to  compute  the  corrective  maneuvers  uses  the  transition  matrix  & . 
Mathematically  the  variation  of  the  satellite  state  can  be  expressed  as: 

^>xx(tf>to)  <&xv(tf>to) 

(7) 

<Pvx(tf,t0)  0vv(tf,to)_ 

The  final  position  of  the  controlled  satellite  using  the  transition  matrix  can  be  formulated  as: 

rf  =  rfref  +  Srf  =  ffref  +  <Pxx(t f  ,t0 )(r0- f0ref  )  +  <PXV (t f  ,t0 )(v0  +  Av0 -  v0ref  )  (8) 

Imposing  the  constraint  that  /y  =  ^fref  ,  the  initial  the  Av0  can  be  obtained  as: 

^^0  ~  ~  ^ xv  ^ xx  0*0  ~  r0,ref  )  ~  ^0  ^ 0,ref  (^) 

In  a  similar  way,  the  velocity  can  be  expressed  as: 

Vf  =  Vfref  +  0vx(tf,to  )(r0  -  r0,ref  )  +  ®vv(tf>  to  )(Vo  +  A^0  ~  V( O.ref  )  (10) 

With  the  constrain  that  v y  +  Av1  =  v ,  the  final  impulse  Av1  is  obtained  by  the  following  equation: 

Av}  =  -  0VX  (r0  -  r0ref  )  -  <PVV  (v0  +  Av0  -  v0ref  )  (11) 
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4.2  Circular  maneuvers  guidance  and  control 

So  far,  all  of  the  RvD  maneuvers  have  been  performed  in  circular  or  near  circular  orbits  due  to  the  fuel 
consumption  and  maneuver  simplicity.  The  formulation  of  these  maneuvers  is  generally  derived  from  the 
Hill-Equations,  which  are  the  simplified  version  of  the  Clohessy-Wiltshire  equation  for  small 
eccentricities  and  relative  distances,  as  shown  in  Eq.  (12). 

2 

x-3co  x  +  2  coy  +  ax 

y  =  -2  cox  +  ay  (12) 

2 

z  -  -co  z  +  az 

In  case  of  circular  orbits,  the  LAREDO  tool  orbital  maneuvers  are  all  based  on  the  Clohessy-Wiltshire 
equations4,  where  the  set  of  maneuvers  implemented  are: 

•  Impulsive  hohmann  transfer:  The  aim  of  this  maneuver  is  to  modify  the  relative  height  of  the  two 
satellites  using  two  impulses  in  the  direction  of  the  Target  velocity.  The  maneuver  duration  is 
equal  to  half  an  orbital  period  (see  Fig.  3). 


Tr^paor*  ui  CBF<he*flH  *i.  XF 


Figure  3.  Impulse  Hohmann  transfer,  Impulse  Hohmann  Rephasing, 
Impulse  Hopping  on  V-bar  and  Continuous  Hopping  on  V-bar  examples. 


•  Impulsive  hohmann  rephasing:  This  maneuver  modifies  the  relative  along  track  distance  of  the 
two  satellites  using  two  impulses  in  the  direction  of  the  Target  velocity.  The  difference  with  the 
Impulsive  Hopping  on  V-bar  is  that  the  maneuver  duration  is  a  multiple  of  the  orbital  period  (see 
Fig.  3). 

•  Impulsive  hopping  on  V-bar:  Maneuver  to  modify  the  relative  along  track  distance  of  the  two 
satellites  using  two  impulses  in  the  direction  of  the  Target  velocity.  The  maneuver  duration  is 
equal  to  half-orbital  period  (see  Fig. 3). 

•  Continuous  hopping  on  V-bar:  The  aim  of  this  maneuver  is  to  modify  the  relative  along  track  shift 
between  the  two  satellites  using  a  constant  acceleration  in  the  direction  normal  to  the  Target 
velocity.  The  maneuver  duration  is  equal  to  an  integer  multiple  of  the  orbital  period  (see  Fig.  3). 

•  Continuous  tangential  transfer:  The  aim  of  this  maneuver  is  to  modify  the  relative  height  of  the 
two  satellites  using  a  constant  acceleration  in  the  direction  of  the  Target  velocity.  The  maneuver 
duration  is  equal  to  an  integer  multiple  of  the  orbital  period  (see  Fig.  4). 
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Figure  4.  Continuous  Tangential  Transfer,  Impulsive  Two  Points  Transfer, 

Station  Keeping  and  Forced  Motion  examples. 

•  Impulsive  two  point  transfer:  The  aim  of  this  maneuver  is  to  take  the  Chaser  satellite  from  an 
initial  point  to  a  generic  final  point,  with  generic  final  velocity  and  with  generic  maneuver 
duration  (see  Fig.  4). 

•  Station  keeping.  The  aim  of  this  maneuver  is  to  keep  the  relative  position  between  the  Target  and 
the  Chaser  as  close  as  possible  to  a  static  reference  position  in  CRF,  using  a  LQ  algorithm.  The 
user  specifies  the  static  reference  position  and  the  duration  of  this  maneuver  (see  Fig.  4). 

•  Forced  motion.  The  aim  of  this  maneuver  is  to  keep  the  relative  position  between  the  Target  and 
the  Chaser  as  close  as  possible  to  a  dynamic  reference  position  in  CRF,  using  a  LQ  algorithm. 
This  dynamic  reference  position  moves  at  a  constant  velocity  in  CRF  from  the  current  position  of 
the  Chaser  to  a  user  defined  final  position.  The  duration  of  this  maneuver,  which  defines  the 
constant  velocity  of  the  dynamic  reference  position,  is  also  specified  by  the  user.  The  trajectory  of 
the  reference  position  is  a  straight  line  in  CRF  (see  Fig.  4). 

4.3  Elliptical  maneuvers  guidance  and  control 

The  Clohessy-Wiltshire  equations  described  in  the  above  section  cannot  be  applied  when  the  orbits  have  a 
significant  eccentricity.  For  these  orbits  a  different  mathematical  approach,  based  on  a  transition  matrix,  is 
applied. 

The  formulation  of  the  elliptical  maneuvers  provided  below  is  described  in  terms  of  time,  but  the 
conversion  to  distance  is  straightforward.  The  elliptical  maneuvers  will  place  the  Chaser  spacecraft  in  the 
same  orbit  as  the  Target  one  but  with  a  certain  time  delay,  tdeiay  ,  with  respect  to  it.  This  means  that  if  the 

maneuver  ends  at  a  time  tnthe  following  relation  must  be  satisfied: 

&f  =xt(tn)-x,  (tn  -  tdelay  )=xc(tn)-  xt  (tn  )  (13) 

The  relative  position  of  the  Chaser  with  respect  to  the  Target  at  the  beginning  of  the  maneuver  can  be 
expressed  in  the  following  way: 


& 0  =xc(to)-xt  (*fl) 


(14) 
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The  Chaser  position  can  be  formulated  as  a  linearization  with  respect  to  the  Target  spacecraft  as  follows: 


&(tn)  =  ®(tn>t0)^0  + 


0  n  0  i 

^ xv,(J  ^ xv,l 


0 


0  n  0  1  0 

_w,U  ^ vv,l  ^ 


xv,n-l 

vv,n-l 


Av0 

Avj 


Av 


n-1 . 


(15) 


Please  notice  that  Avn  is  omitted  since  the  position  of  the  satellite  at  the  end  of  the  set  of  maneuvers  is  not 

modified.  Using  the  Lagrange  operators  X  ,  the  function  to  be  minimized  to  perform  the  maneuvers  in  an 
optimal  way  is: 


/  =  Trace(  Av  T  Av  )  +  Trace(  X 


-<p(tn,t0)&0 


0  n  0  i  0  i 

^ xv, (J  ^ xv, 1  ^  xv, n-1 

0  n  0  i  0  i 

^ vv,(J  ^ vv,l  **•  ^ vv,n-l 


Av 


)m 


Minimizing  the  last  equation  with  respect  to  Av  and  A  : 


%= 

dA 


df  -  o  -  2Av 

0xv,O 

0  , 

^  XV, 1 

0  /I 

^  xv, n-1 

d(Av ) 

0W,O 

0vv,l 

0  1 

^vv,n-l  _ 

0  =  Scf  -0(tn,to)Sco  - 

>rv,0 

0w,O 

0  i  0 

0  i  0 

^ vv,l  •" 

A 


xv, n-1 
vv,n-l 


Av 


(17) 


(18) 


Putting  these  last  equations  in  a  matrix  form,  and  solving  to  obtain  the  multiple  “n”  Av  necessary  to 
achieve  the  required  position  at  time  tm  the  following  equation  arises: 


~2I3x3 

°3x3 

°3x3 

~0xv,O 

1 

1 

-1 

Av~ 

°3x3 

21 3x3 

°3x3 

-  0 XV, 1 

-0  i 

^ vv,l 

°3nxl 

A  _ 

°3x3 

°3x3 

-  2I3x3 

-  0 xv, n-1 

~0vv,n-l 

_&f  ~  0(tn<tO  )<&o 

0xv,O 

0xv,l 

-  0 xv, n-1 

°3x3 

°3x3 

1 

0vv,l 

-  0vv,n-l 

°3x3 

°3x3 

4.4  Orbital  navigation 

The  orbital  navigation  is  carried  out  using  measurements  from  ground  station,  but  also  with  the  modeled 
satellite  sensors,  accelerometers,  IMU,  RF,  camera,  and  lidar. 

The  navigation  solution  is  obtained  using  a  Kalman  filter  formulation. 
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5.0  ATTITUDE  MANOEUVRES 

The  Attitude  plays  a  very  important  role  when  the  final  part  of  the  RvD,  the  docking,  is  carried  out.  The 
attitude  needs  to  be  accurate  enough  to  allow  the  satellites  to  dock.  The  spacecraft  must  be  able  to  navigate 
either  when  only  inertial  sensors  (e.g.  Star  tracker)  are  available  or  when  relative  sensors  are  available. 
The  guidance  in  these  two  cases  is  completely  different. 

5.1  Attitude  maneuver  guidance 

The  satellite  implements  a  few  sets  of  guidance  profiles  to  perform  the  required  inertial  and  relative 
maneuvers.  The  type  of  guidance  will  depend  upon  the  available  sensors  and  upon  the  satellite  relative 
distances.  Within  LAREDO,  the  following  guidance  profiles  are  considered: 


Figure  5.  Fixed  Relative. 


Figure  6.  Target  Pointing. 


Figure  9.  GS  pointing. 


•  Fixed  relative  attitude.  Guidance  to  maintain  the  Chaser  in  a  constant  orientation  with  respect  to 
the  Target  (see  Fig.  5  ). 

•  Target  pointing  (see  Fig.  6).  Guidance  to  point  towards  the  Target  satellite.  Very  useful  if  the 
Chaser  satellite  is  desired  to  point  always  towards  a  sensor  placed  on  the  Target  spacecraft  (e.g. 
lasers  or  cameras). 

•  Sun  pointing.  Guidance  to  point  towards  the  Sun.  Useful  mode  if  batteries  are  desired  to  be 
charged  or  if  Collision  avoidance  maneuver  is  activated  (see  Fig.  7). 

•  Nadir  pointing.  Guidance  to  point  towards  the  Earth  Center  (see  Fig  8). 

•  GS  pointing.  Guidance  to  point  towards  the  Ground  Station  (GS).  This  type  of  guidance  is 
required  if  the  antenna  is  not  directional  when  uploading  or  downloading  information  from  GS 
(see  Fig.  9). 

5.2  Control 

The  control  used  in  attitude  is  a  Proportional  Derivative  Controller  (PD) 5.  This  type  of  attitude  control  is 
simple,  but  its  performance  is  very  good.  The  equation  of  the  implemented  controller  is  as  follows: 
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TCx  ~2Kx  cllE(l4E  +  Kxd  wEx 

Tcy  -  2Ky  q2E(i4E  +  Kyd  wEy  (20) 

TCz  =  2Kz  cl3E(l4E  +  Kzd  wEz 

Where  Kx,  Ky  ,KZ  are  the  proportional  gains  of  the  x,  y,  z  axis;  Kxd,Kyd,Kzd  are  the  derivative 
gains  of  the  x,  y,  z  axis;  q1E  ,  q2E  and  q3E  are  the  vectorial  part  of  the  quaternion  error  (error  between 
the  desired  and  actual  attitude);  q4E  is  the  scalar  part  of  the  quaternion  error;  wE  ,  wE  and  wE  are  the 
x,  y  and  z  components  of  the  angular  velocity  error. 

5.3  Attitude  navigation 

Attitude  navigation  is  carried  out  autonomously  since  no  measurements  from  the  ground  are  required.  The 
sensors  modeled  are  Star  trackers,  RF  system,  Camera  system,  Gyroscopes  and  Sun  sensors.  The  sensor 
noises  are  filtered  using  a  Kalman  Filter6  to  obtain  a  clear  signal. 

6.0  TEST  CASES 

This  paper  reports  the  simulation  results  obtained  for  a  Mars  Sample  Return  mission  scenario7.  The 
simulation  of  the  this  mission  scenario  includes: 

•  Launch  ascent  trajectory  optimization  of  the  PAV  to  insert  it  in  a  circular  orbit. 

•  A  change  of  plane  maneuver  to  place  both  satellites  in  the  same  orbital  plane. 

•  RvD  phase  of  the  satellites  in  circular  orbits. 

Further  more,  the  paper  also  shows  the  results  of  a  RvD  in  elliptical  orbits  to  illustrate  the  tool’s 
performances.  The  attitude  guidance  and  control  is  already  presented  in  the  Attitude  Guidance  Maneuvers 
section  (Fig.  5,  Fig.  6,  Fig.  7,  Fig  8  and  Fig.  9).  These  figures  are  obtained  directly  from  the  LAREDO, 
where  all  the  sensors,  and  actuators  were  activated. 

6.1  Mars  sample  return  launch  ascent  results 

The  launch  is  performed  from  45  degrees  latitude  and  0  degrees  longitude  on  the  Mars  surface.  In  this 
simulation,  the  atmospheric  and  non-sphericity  perturbations  are  taken  into  account.  The  launch  ascent 
vehicle  characteristics  and  the  desired  final  orbit  is  provided  in  the  tables  shown  bellow  (Table  1  and 
Table  2).  Note  that  the  Target  spacecraft  is  also  simulated,  and  its  initial  argument  of  latitude  is  78.70  deg. 

Table  1.  Launch  ascent  Final  orbit 


Semi-major  axis  (km) 

3897.514 

Eccentricity 

0 

Inclination  (deg.) 

45 

Right  ascention  node  (deg) 

265.811 

Argument  of  perigee  (deg) 

253.629 

Max.  relative  Argument  of  latitude  (deg) 

4 

Min.  relative  Argument  of  latitude  (deg) 

2 
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Table  2  Launch  ascent  Rocket  characteristics 


First  Stage 

Coast  Arc 

Third  Stage 

Reference  diameter  (m) 

0.563 

0.563 

0.563 

Nominal  burn  time  (s) 

197.566 

- 

672.197 

Propellant  mass  (kg) 

360 

- 

116 

Gross  mass  (kg) 

456.6 

- 

213.8 

Specific  Impulse  (s) 

306 

- 

325 

Fig.  10  and  Fig.  11  show  that  the  circular  orbit  of  500  km  height  is  achieved.  The  atmospheric  pressure 
losses  can  be  clearly  observed  in  the  thrust  profile  of  Fig.  12,  where  the  thrust  increases  when  higher 
altitudes  are  reached.  The  rocket  loss  of  mass  is  represented  in  the  Fig.  13.  This  figure  clearly  shows  that 
the  mass  is  reduced  when  the  first  stage  of  the  rocket  is  deployed.  With  the  optimized  trajectory,  the  rocket 
can  deliver  10.32  kg  of  payload  mass,  deploying  it  at  a  relative  argument  of  latitude  with  respect  to  the 
Target  of  2  degrees. 


Figure  10.  Launch  ascent  height. 


0  1W3  200  500  *00  6W  TOW  6O0  9CW 
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Figure  11.  Launch  ascent 
apocentre  and  pericentre. 


Figure  12.  Launch  ascent  Thrust. 


Figure  13.  Launch  ascent  Mass. 


6.2  Mars  sample  return  change  of  plane 

Generally,  the  launcher  does  not  deploy  the  PAV  exactly  in  the  desired  orbit.  The  PAV  is  delivered  into  a 
circular  orbit  that  is  two  degrees  different  in  inclination  and  Right  ascension  with  respect  the  Target  Orbit, 
and  20  km  below  (see  Table  3).  Therefore,  a  change  of  plane  maneuver  is  performed. 
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Table  3.  Change  of  Plane  orbital  elements 


Initial  orbit  before  change  of 
plane 

Final  orbit  after  change  of 
plane 

Semi-major  axis  (km) 

3877.515 

3897.514 

Eccentricity 

0 

0 

Inclination  (deg) 

45 

47 

Right  ascension  node  (deg) 

265.811 

268.811 

Argument  of  perigee  (deg) 

253.629 

253.629 

True  anomaly  (deg) 

121.035 

123.035 

Kid6, 


K  m 


Figure  14.  Change  of  Plane  Maneuver. 

6.3  Mars  sample  return  RvD  in  circular  orbit 

Once  the  change  of  plane  is  performed,  the  RvD  can  be  done  nicely  in  the  coplanar  circular  orbits.  The 
satellites  starts  from  the  change  of  plane  maneuver  at  a  relative  height  of  20  km  and  -60  km  of  along 
track,  performing  the  required  maneuvers  (also  applying  corrective  maneuvers)  until  the  docking  is 
achieved  (see  Fig.  15  and  Fig.  16).  The  maneuver  Timeline  is  as  follows: 

•  Impulsive  Hohmann  Transfer  to  place  the  satellites  at  the  same  height 

•  Impulsive  two  point  transfer  to  place  the  Chaser  -2.7  km  in  along  track 

•  Station  keeping  at  -2.7  km  along  track  for  1000  seconds 

•  Impulsive  Hopping  on  V-bar  to  place  the  Chaser  -300  km  in  along  track 

•  Station  Keeping  at  -300  km  along  track  for  1000  seconds 

•  Forced  Motion  for  2000  seconds  to  place  the  Chaser  -100  km 

•  Forced  Motion  for  2000  seconds  until  final  Docking. 
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Figure  15.  Rendezvous  and 
Docking  trajectory. 


Figure  16.  Rendezvous  and 
Docking  position. 


Figure  17.  Rendezvous  AV. 


6.4  RvD  in  elliptical  orbit 

This  test  case  is  presented  to  show  the  results  obtained  using  the  RvD  algorithms  implemented  for 
elliptical  orbits.  The  initial  satellite  states  used  for  this  simulation  are  provided  in  the  table  shown  below 
(see  Table  4): 


Table  4.  Elliptical  RvD  orbital  elements 


Chaser  Orbital  Elements 

Target  Orbital  Elements 

Semi-major  axis  (km) 

9257.165 

9306.165 

Eccentricity 

0.607 

0.607 

Inclination  (deg) 

45 

45 

Right  ascension  node  (deg) 

232.934 

232.934 

Argument  of  perigee  (deg) 

345.076 

345.076 

True  anomaly  (deg) 

-0.233 

-0.033 

The  initial  relative  position  of  the  Chaser  with  respect  to  the  Target  spacecraft  is  a  difference  of  20  km  in 
height  and  nearly  13  km  in  along  track.  It  can  be  observed  from  the  table  above  (Table  4)  that  both  the 
Chaser  and  Target  spacecrafts  are  in  high  elliptical  orbits.  The  Elliptical  RvD  maneuvers  Timeline  is  as 
follow  (see  Fig.  18): 

•  Elliptical  Transfer  to  place  the  satellites  at  the  same  height 
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•  Elliptical  Re-phasing  in  Time.  Chaser  is  placed  at  1 -second  delay  with  respect  to  the  Target. 

•  Elliptical  Hopping  on  V-bar  to  place  the  Chaser  -  5km  in  along  track. 

•  Elliptical  Station  Keeping  at  -5km  along  track. 

•  Elliptical  two  points  transfer  to  place  the  Chaser  -1  km  in  along  track. 

•  Elliptical  Forced  Motion  until  Docking  is  performed. 


EKdcal  Rccrftavng  in  Tm. 
EK*>caJ  Mopping  on  V-bar 
EM<*SKwng<V-b*(l 
EKfecai  Two  Point  Transltr 
Elfcal  F  Mobon/S  Katpng 


’’  V  2  4  6  8  1*0  12 

Ttmew  Kl<J* 


Figure  18.  Elliptical  Rendezvous. 


7.0  CONCLUSION 

LAREDO  is  a  very  sophisticated  tool  to  perform  Launch  ascent  trajectory  optimizations  and  RvD  for 
elliptical  and  circular  orbits.  The  launch  ascent  optimizer  provides  a  very  accurate  way  to  obtain  an 
optimized  trajectory  when  launching  from  the  Moon,  Mars  and  the  Earth.  For  the  RvD  phase,  Venus, 
Deimos,  Europa  and  Triton  can  also  be  used  as  central  bodies  in  addition  to  the  Earth,  the  Moon  and  Mars. 
The  RvD  in  circular  orbits  uses  the  Clohessy-Wiltshire  equations  with  the  aid  of  corrective  maneuvers. 
The  RvD  in  elliptical  orbits  is  based  on  the  transition  matrix  to  cope  with  the  non-linearity  of  the  problem. 

The  results  presented  in  this  document  show  the  good  performances  of  the  launch  ascent  trajectory  for  a 
two-stage  rocket.  The  simulations  presented  for  RvD  are  performed  with  all  perturbations  acting  upon  the 
satellite  and  using  sensors  and  actuators.  Both  the  Orbital  and  Attitude  algorithms  implemented  in 
LAREDO  are  extensively  proved  and  documented. 

Future  work  will  be  done  to  prove  the  orbital  and  attitude  algorithms  in  a  test  bench  to  see  their 
performances  in  a  more  realistic  way. 
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Laredo  Description  (1) 


what's  Laredo? 

o  SW  tool  designed  and  developed  by  GMV  S.A.,  Spain,  under 
GSTP  contract  with  ESA/ESTEC. 

o  Devoted  for  mission  analysis  and  design  activities  in  order  to 
define/verify  the  system  requirements  (from  orbit/attitude 
control  point  of  view)  of  a  dedicated  RdV/FF  mission,  helping 
on  the  preparation  of  future  control  and  operations  required 
capabilities. 


www.gmv.com 


Laredo  Description  (1) 


What  Laredo  is  meant  for? 

o  Cover  launching  phase  of  the  spacecrafts  involved  in  the 
RvD/FF  once  it  has  left  from  Earth.  Moon  and  Mars  surfaces. 

o  Simulate  Orbital  Rendezvous  and  Dockinq/Formation 
Flying  sequences  on  Circular  and  Elliptical  orbits  for 
exploration  mission-scenarios  such  as:  Earth.  Venus.  Moon. 
Mars.  Deimos.  Triton  and  Europa. 

o  Simulate  Orbital  Rendezvous  and  Docking  sequences  on  co- 
linear  libration  points  in:  Sun  +  EarthXMoon  Barycentre.  Earth 
+  Moon  and  Mars  systems. 


o  In  particular  Earth  Observation  missions  managing  several 
satellites  under  a  particular  configuration  (Rdv/FF)  is  a  target 
scenario  valid  in  Laredo. 
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Laredo  Description  (2) 


Laredo  first  level  apabilities  (1/2): 

o  Open  loop:  no  sensors  or  actuators 
o  Corrective  manoeuvres 
o  Three  degrees  of  freedom 

o  Environment:  drag,  solar  radiation,  third  body,  non-sphericity 
o  Launch  ascent  optimisation 
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Laredo  capabilities  (2/2): 


o  Available  Orbital  Manoeuvres  for  Rendezvous/FF: 


Circular: 

Elliptical: 

>  Hopping  on  V-bar 

>  Hopping  on  V-bar 

>  Re-phasing  in  Time 

>  Hohmann  Re-phasing 

>  Height  Transfer 

>  Two  Points  Transfer 

>  Hohmann  Transfer 

>  Applied  Delta-V 

>  Two  Points  Transfer 

>  Station  Keeping 

>  Forced  Motion 

>  Continuous  Tangential  Transfer 

Libration  points: 

>  Continuous  Hopping  on  V-bar 

>  Hopping  on  V-bar 

>  Height  Transfer 

>  Applied  Delta-V 

>  Two  Points  Transfer 

>  Station  Keeping 

>  Applied  Delta-V 

>  Station  Keeping 

>  Forced  Motion 

>  Forced  Motion 
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Laredo  Description  (4) 

Laredo  second  level  capabilities: 

o  Close  loop:  feedback  of  sensors  and  actuators 
o  Six  degrees  of  freedom  (translation  +  attitude) 
o  Communications  module 
o  Guidance,  Navigation  and  Control  (GNC) 
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Introduction  to  Laredo  Execution 
o  Definition  of  inputs 
o  Execution  of  simulation 
o  Visualization  of  results 
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General  Architecture  (2) 


The  architecture  of  Laredo  has  been  design  to 
achieve  these  three  tasks. 
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)  FAM0SV2  vl.5 


Input  Definition  ■  Simulation  Run  Outputs  Help 


General  Architecture  (3) 


^jnjxj 


Simulation  parameters 

|  atv.dat  ~^\ 

S/C  Layout 

|  atv.dat  ~^\ 

RvD  Initial  States 


|  atv.dat 

Timelines 

|  atv.dat  3] 

Launching  options 

|mars_arg.dat 

Optimization  param. 

|mars_arg.dat 

Rocket  parameters 

|  mars_mass.dat 

PAV/Or biter  Orbits 

|  mars_arg.dat  ~^\ 

Change  of  Plane 

|  mars_cop.dat  ~^\ 

”3 


LQR  gains 

|  lqr.dat 


Inputs  Definition 

atv.dat 

Simulation  Parameters 

Step  size  |  IQ 
Simulation  time  |  21000” 


~3 


|  Modified  Julian  Day  2000 

Initial  epoch  |  g~~ 


Central  body  |  Earth 


Aero  drag  |  Computed 

Atm  selection:  |MarsGram1  SS0 
MCD  Folder:  |" 


VEMCD_DATA/ 


Solar  radiation  |  Computed 


Non-sphericity  |  Computed 

Number  of  zonals  |~" 


20 


1  0 

Number  of  tesserals 

Third  Body  Perturbation 

20 

|  Ignored 

d 

Uranus  |gn0red 

zl 

Nereid  ||gnored 

d 

|  Ignored 

Id 

Neptune  |gn0red 

d 

Deimos  |  |gn0red 

d 

(ignored 

d 

|  Ignored 

d 

phobos  ||gnored 

d 

|  Ignored 

d 

Moon  Computed 

d 

1°  |  Ignored 

d 

|  Ignored 

d 

^un  |  Computed 

d 

Ganymede  ||gnored 

d 

|  Ignored 

d 

Callisto  |  Ignored 

d 

Input  Data 
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MY 


General  Architecture  (4) 


+)  FAM0SV2  vl.5 


Input  Definition  Simulation  Run  ■  ©otpoti  .  Helh 


^jnjxj 


Vi  .5 


Simulation  Run 

i 

atv.dat 

Run! 

Abort 

Project  Definition 

Projects 

|  atv.dat 


~3 


Launch  Window 

|  nnars_mass.dat 


Project  definition:  |  Rendezvous  Manoeuvres 
Simulation  Parameters  file: 


atv 


'if 


atv 


atv 


Initial  guess  from:  |  Previous  output 
□  rbiter  Layout  file: 

Launching  Options  file: 

Optimization  Parameters  file: 

Rocket  Parameters  file: 

PAV/O  rbiter  Orbits  file:  I 


~3 


atv 


atv 


atv 


■=1 
'd 
d 

-J  Change  of  Plane  file  |"~ 

'd 


Timeline  file: 
LQR  gains: 
Change  of  Plane? 


atv 


Iqr 


No 


~3 


'd 

'd 


Run 

Simulation 


PAV  and  0  rbiter  roles 
(for  RvD  only) : 


Ignore  launcher  optimization:  Input  from  initial  states  file 


3 
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MY 


?  FAM0SV2  vl.5 


'  InpnttJefinition  Simulation  Run  Outputs  Help 


Rend,  fc  Dock. 


Available  Projects 

last 


iJ 


Delete  folder 


RvD  Dynamics 

[Trajectory’  X-Z 

Lissajous  Dynamics 

j  Position  YX 

Budgets 

jAccel.  in  C RF  ^~| 

Output  text  files 

j  Manoeuvres  list 


-1000  r 


0  - 


1000  - 


2000 


3000 


4000  - 


5000  - 


6000 


7000 


0.5 


General  Architecture  (5) 


-Ini 

J 

a - 

vl  .5 

Trajectory  in  CRF  (height  vs.  X) 


- i 

l - , - 

- - 1 - 1 

l - l 

l - 1 

- 1 

_ rfi _ 

_ j 

i_ 

j 

_ j 

_ _ j 

r 

r 

_ 

r\7 

-  Free  Drift 

-  Impulsive  Hohmann  Transfer 

O  Station  Keeping  (user  defined) 

-  Impulsive  Hopping  on  V-bar 

O  Station  Keeping  (user  defined) 

-  Forced  Motion 

O  Station  Keeping  (user  defined) 

-  Forced  Motion 

-  Free  Drift 

- 1 - 

“  Y  “  1 

- i - 

- 1 - 

- 1 

— 

- 1 

_ j _ 

j 

_ _ j 

_L 


-0.5 


-1.5  -2 

X  (m) 


-2.5 


-3.5 


x  10 


FOMS 
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Laredo  Propagator  (1) 


o  Full  Cowell  numerical  propagator, 
Domand-Prince  (ode5) 

•  Propagation  around  the  following 
central  bodies: 

>  Earth,  Venus,  Mars,  Moon,  Deimos, 

Europa  and  Triton 

•  Four  types  of  perturbation  can  be  considered: 

>  Aerodynamic  drag  perturbation 

>  Solar  radiation  pressure  perturbation 

>  Non-sphericity  perturbation 

>  Third  body  perturbation 

•  Simplified  propagator  for  launcher  ascent: 

>  J2  and  Aerodynamic  drag  perturbations 


•  Attitude  propagation 


>  ^och  MJD2000 

S  C_Posrtion_MEE2DDD 

>  S  C_Atitude  _MEE2000 

>  Control_Acc_MEE2DDD 

SC_\^locity_MEE2DDD 

>  Planets  Positions 

>  Sun  Position 

Ajx_V£rs 

The  Propagator 
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o  Atmosphere  modelling 


•  Earth 

>  US  Standard  Atmosphere  1976  tabulated  data  below  20 
km. 

>  GRAM99  model  from  27  to  90  km. 

>  MSIS-86  model  above  120  km. 

>  Fairing  between  the  appropriate  models  from  20  to  27 
km  and  from  90  to  120  km 

•  Venus 

>  Pioneer  Venus  Thermospheric  Model  (1983) 

•  Mars 

>  MARS  GRAM  atmosphere  model  developed  by  NASA  in 
1989  or  MCD  (Mars  Climate  Database) 

•  Triton 

>  Exponential  approach  has  been  considered  for  the 


density  calculation. 
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o  Solar  Radiation  Pressure  Perturbation 

•  Cannon-ball  model 

•  Eclipse  by  central  body  considered 


o  Non-sphericity  Perturbation 

•  Venus:  2x2  (NASA  technical  memo) 

•  Earth:  50  x  50  (GEMT3) 

•  Moon:  4x4  (Ferrari  1977) 

•  Mars:  50  x  50  (Goddard  Mars  Model  -  2B) 


o  Third  Body  Perturbation 

•  Sun,  planets  and  major  satellites  of  available  central  bodies. 

o  Solar  System  Bodies  Ephemerides 

•  Sun,  planets  and  the  moon:  JPL  planetary  and  lunar 
ephemerides 

•  Natural  Satellites  (except  Moon):  Analytical  expressions  for 
natural  satellites  ephemerides  (1987) 
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Sensors  &  Actuators 


o  Sensors 

•  RF 

•  CAMERA 

•  LIDAR 

•  IMU 

•  STR 

•  SAS 

o  Actuators 


MEEHIECRF 

CD- 


C  Poflfoi  CRF  rea 


C  AM  q  I.IEEUE  real 

CD— 

CRF2MEEZH 

-+ 


A_mal 


Galt  Imc?  leclcr 


roL^al 


H_mal 


AB_mal 


R0F2I.I  EEZHI 


HEEZm-CRDF 


V_a 

VJ 

'j_aj 


REl  P05 C  P.B  F 


(D  D  D> 


C  COM  CR0F 


Hr.lkfirjM_citrf 


□ 


■«>» 

SC2LIDAR 


t  -+o 


FlIFtf  LIPAP. 


9ccpe2 


RElPof.Udajir 
RflP«_Uda 

RElPoi.Udsjull 

Error  flddldoi 


LIDAR2SC 


Hr.lffir_p«_citJf 


T  m  jnal  I  ■ 

tUfc|+ 

- ►ajk _ I  , 


(D  n  i) 


c«(Ur.t«i) 


fblufjdr 


XD 


rag 


i 

Eiiife 


Re&foial 

Operator 


-KD 


C  PoflUdk  CRF  teifor 


diflt 


DertettK 


C_Ve  bctV_C  R  F_se  isor 


Attitude  Thrusters 
Orbital  Thrusters 
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Fix  Relative 


Target 

Pointing 


Attitude  Guidance 
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o  Navigation 

•  Orbital  and  Attitude  Estimation  is  performed  with  Kalman  filter. 

•  Filtering  of  different  sensor  measurements. 

•  Propagation  of  state  is  performed  when  no  measurements  are 
available. 

o  Communications 

•  No  limit  in  GS  number 

•  GS  visibility  computed. 

•  Signal  delay  simulated 
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o  Circular: 


Validation  Scenarios 


Mars  Sample  Return  (Launch  ascent  from  Mars) 

Mars  Sample  Return  (Change  of  Plane  &  Rendezvous) 


Deimos  Sample  Return  (Rendezvous) 
Human  Lunar  Mission  (Ascent) 

Human  Lunar  Mission  (Rendezvous) 
ATV  Mission  (Rendezvous) 

Combined  Earth  Observation  (FF) 


S^uiKh^I] 


Courtesy  of  NASA 


o  Elliptic 

•  Launch  Ascent  to  Elliptical  Orbit 

•  ATV  RvD  with  Corrective  Manoeuvres 

•  Elliptical  Approach 


Courtesy  of  ESA 


www.gmv.com 


